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Abstract. The simultaneous beta decay of two neutrons in a nucleus without the emission of 
neutrinos (called neutrinoless double beta decay) is a lepton number violating process which is 
not allowed in the Standard Model of particle physics. More than a dozen experiments using 
different candidate isotopes and a variety of detection techniques are searching for this decay. 
Some (EXO-200, Kamland-Zen, GERDA) started to take data recently. EXO and Kamland- 
Zen have reported first limits of the half life T®J 2 for 136 Xe. After a decade of little progress 
in this field, many new results will soon scrutinize the claim from part of the Heidelberg- 
Moscow collaboration to have observed this decay. The sensitivities of the different proposals 
are reviewed. 



1. Introduction 

For 35 isotopes, /3 decay is energetically forbidden but double beta decay (2uf3f3) is allowed 



This process has been observed directly for 11 isotopes with half lives between 7 • 10 yr and 



Since neutrinos have no electric charge, there is no known symmetry which forbids that they 
mix with their anti-particles. As a consequence double beta decay without neutrino emission 
(Oz^/3/3) might occur as well which is predicted by almost all extensions of the Standard Model 
of particle physics (see e.g. [HE]). Its observation would imply that lepton number is violated 
since only electrons are emitted. 

The experimental signature of Ov(3(3 is a line at the Qpp value of the decay if the sum of the 
electron energies is histogrammed. 

Part of the Heidelberg- Moscow collaboration claims to have evidence for this line for 76 Ge 
with Ty 2 = (1.191q23) ' y r [6j at a 4<r level. In a more sophisticated analysis a significance 
beyond 6a was calculated [7J. The current experiments using 136 Xe (EXO-200 and Kamland- 
Zen) and 76 Ge (GERDA) will scrutinize this result within the next 12 months. 

In total more than a dozen large scale experimental programs are suggested or under 
construction to search for Ov/3/3. These programs are compared in this article and also the 
status of theoretical matrix element calculations is discussed. The latter are needed to convert 
an experimental measurement or limit on Ty 2 to a particle physics parameter. For general 
reviews the reader is referred to the literature [8j SJ EJ [5] . 

There are also other related processes like double positron decay or double electron capture. 
While 0v/3/3 is already a suppressed process, the other decays are expected to be even rarer 
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unless there is some resonance enhancement [TUJ, [TTJ [TJ1 El HB]- In this article only 0^/3/3 decay 
searches are discussed. 



2. Motivation for Oi//3/3 

Since neutrinos have mass but no electric charge, there is no known symmetry which forbids 
additional terms in the effective Lagrangian besides the standard Dirac mass term mo 



L Yu k = rn D v L v R + -m L v L {v L ) c + -m R {v R ) c v R + h.c. (2) 
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The subscript L stands for the left-handed chiral field ul = ^(1 — 75) v and R for the right- 
handed projection |(l + 7s) v. The superscript C denotes charge conjugation, i.e. v c stands for 
an (incoming) anti- neutrino. The m R term describes therefore an incoming neutrino v R and an 
outgoing anti-neutrino (i / R ) c , i.e. this term violates lepton number by 2 units. The eigen states 
of the mass matrix are of the form {y + v c ). Consequently neutrinos are expected to be - in 
general - their own anti-particles, i.e. Majorana particles. 

Neutrinos (or anti-neutrinos) are produced in charged weak current reactions and - depending 
on the charge of the associated lepton - only one chiral projection couples ("V-A" current of 
weak interactions). For example in /3 decay n — > pe~ v e ,R, a right-handed anti-neutrino couples 
which can be decomposed in the mass eigen states V{ and helicity {h = op) eigen states: 

1 3 • 

Ve,R = V e -{1 + 75) = X U ei(Vi,h=+l + ~Vi,h=-l) (4) 
- i=l 

Here, U is the PMNS mixing matrix [15], rrii are the mass eigen values, and E is the neutrino 
energy. 

For Dirac particles, only detection reactions like p v e ^ R — > n e + are possible. If, on the other 
hand, neutrinos are massive Majorana particles, then the helicity suppressed component Uih=—i 
of v e ,R can undergo the reaction nv e ^ — >• pe~ . Here 
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Ve,L = ^(1 - 75K = V] Uei(Vi.h=-l + '-^ v i,h=+l) (5) 
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Taking both processes involving neutrons together we have 2n — > 2p + 2e~ (or better 
(A, Z) — > (A, Z + 2) + 2e~) mediated by a massive Majorana neutrino with an effective coupling 
strength which is called the Majorana mass: 
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m t3P = \Y. U ei- m i\ ( 6 ) 
i=l 

The helicity suppression (rrii/E) 2 , which is e.g. 10~ 14 for a neutrino mass of 0.1 eV and a 
neutrino energy of 1 MeV, is compensated by the large number of nuclei per mole. 

Schechter and Valle showed that the observation of ensures that neutrinos have a 

Majorana component [16]. Recently it was pointed out that the "guaranteed" Majorana mass 
through radiative corrections is however only in the range of 10 -24 eV [T7], i.e. negligible 
compared to the mass scales of neutrino oscillations [15] . Consequently, if Of/3/3 is observed, 



other mechanisms like the exchange of supersymmetric particles or heavy Majorana neutrinos 
might be the dominating process and the known neutrino could even be (effectively) a Dirac 
particle. While for the initially motivated light neutrino exchange the coupling strength is 
proportional to (Eq. [6]), the one for processes with heavy fermion exchange is proportional 
to 1 /Mp with Mp being the mass of the exchange particle (provided that the helicity of the two 
leptonic currents are the same and M F > q 2 )^\ In such processes, lepton number violation or 
lepton flavor violation can be accessible with accelerator experiments as well [18j . If lower limits 
on Mp from e.g. LHC are higher than typically 10 TeV, then the contributions to Ov/3/3 become 
smaller than the ones expected from light neutrino exchange with a mass of 0.05 eV [3]. The 
argument can be turned around: 0u/3/3 will provide also information on TeV scale physics. 

The exchange of light neutrinos is discussed predominantly. Its strength depends on U which 
can be parameterized by 3 rotation angles and 1 phase (#12, #13, 6*23) and 5, all measurable by 
neutrino oscillation experiments), and 2 additional phases («2i and 031). The latter are called 
Majorana phases and influence processes like Oz^/3/3. Oscillation experiments measure difference 
of squared masses (Am^ = m\ — mf = (7. 581^26) ' 10 -5 eV 2 and |Am3i| 2 = |m 2 — m\\ = 
(2.35±g;J|) • 10~ 3 eV 2 [15]). Knowledge on the absolute mass scale comes from e.g. beta decay, 
Qvf3fl or cosmology. 

One can estimate possible values of tjirr as a function of the lightest neutrino mass given 
the experimental knowledge about Am^ and the measured rotation angles and allowing the 
Majorana phases to take any value. Three cases can be discriminated |19j . 

• If the lightest mass is larger than all Am^, all m; are similar (degenerate masses). 

• If m\ is the smallest mass (normal mass hierarchy) the 3 terms of mM can cancel for 
mi < J Am^i- For m\ ~ typical values are a few meV. 

• If 1JI3 is the smallest mass (inverted hierarchy, Am^ < 0) there exists for mog a lower 
bound of (19l^ 5 ) meV and an upper bound of 50 meV [2]. 

If additional sterile neutrinos exist, i.e. if there are additional terms in Eq. [6j then even for the 
inverted hierarchy no lower bound exists |20l |2~T] . 

In summary, 0v/3(3 - together with other input from e.g. neutrino oscillation and LHC 
experiments - provides an important window to extensions of the Standard Model with lepton 
number violation. 

3. Experimental sensitivity on T®J 2 

All (but one) searches have so far only observed event counts in the region of interest around 
Qpp which are consistent with the expectation from background Abkg- This number - if it scales 
with the detector mass M - is given by 

Abkg = M ■ t ■ B ■ AE (7) 

Here t is the measurement time, B is the so called background index given typically in 
cnts/(keV-kg-yr), and AE is the width of the search window which depends on the experimental 
energy resolution. Note that this equation is only an approximation. Experiments normally 
take the (nonlinear) shape of the background spectrum in a fit into account. 

The non-observation is converted to an upper limit on the number of signal events A a i g which 
is related to the half life Ty 2 of a given isotope A by 

A sig = In 2 • N Avg • e ■ V ' M ■ t/(m A ■ lty 2 ). (8) 

To be more precise: the amplitude has dimension mass' " 5 = Mg 4 • Mp 1 with M B being the mass of a scalar or 
vector exchange particle like the W boson. 



Table 1. List of phase space factor (Eq. 10), Qpp, natural abundance of /3/3 isotope, half life 
of 2^/3/3, and experiments for the most interesting Oz^/3/3 isotopes. Half lives for 2v(3(3 are taken 
from [U Qpp for 136 Xe from |22| . and all other numbers from [3]. The G 0l/ values have been 
scaled to the same nuclear radius and qa coupling. 



isotope 


I 10 " 14 ] 

L yr J 


Qpp 

[keV] 


nat. ab. 

[%] 


rp'iv 
1 l/2 

[10 20 y] 


experiments 


4 «Ca 


6.3 


4273.7 


0.187 


0.44 


CANDLES 


76 Ge 


0.63 


2039.1 


7.8 


15 


GERDA, Majorana Demonstr. 


82g e 


2.7 


2995.5 


9.2 


0.92 


SuperNEMO, Lucifer 


100 Mo 


4.4 


3035.0 


9.6 


0.07 


MOON, AMoRe 


116 Cd 


4.6 


2809.1 


7.6 


0.29 


Cobra 


130 Te 


4.1 


2530.3 


34.5 


9.1 


CUORE 


136 Xe 


4.3 


2457.8 


8.9 


21 


EXO, Next, Kamland-Zen 


150 Nd 


19.2 


3367.3 


5.6 


0.08 


SNO+, DCBA/MTD 



iVAvg is the Avogadro constant, e the signal detection efficiency, r\ the mass fraction of the 0v(3(3 
isotope, and rriA the molar mass of the isotope. 

If Abkg < 1 the experimental sensitivity scales with M ■ t while for Abkg >> 1 the e.g. 90% 
C.L. limit on the half life (assuming there is no signal) is given by 

If systematic errors become important e.g. if the energy resolution or the spectral background 
shape is not well known, then the sensitivity is reduced. 

4. The nuclear matrix element 

(Ty 2 ) is the product of three factors: a phase space factor G 0l/ , a nuclear matrix element M 0u 
and a particle physics factor. 

If the exchange of light Majorana neutrinos is dominating, Ty 2 for a given isotope A is |23] 

2 

[T°J 2 (A)]- 1 = G°"(Qw,Z) ■ \M°»(A)\ 2m f (10) 

m e 

Here, m e is the electron mass. Values for G 0u are listed in Tab. [TJ Obviously, M 0u (A) is 
needed to compare results from different isotopes or to extract information about the particle 
physics parameter. 

The M 0u calculations are difficult and can only be done using approximations. Traditionally, 
interacting shell model (ISM) [241 [25] and quasi particle random phase approximation (QRPA) 
calculations have been performed |27[ I33j . Recently new approaches like the interacting boson 
model (IBM) (29] [301 EI] > the generating coordinate model (GCM, also called energy density 
functional EDF) |28j and the projected Hartree-Fock-Bogoliubov (pHFB) method [32] have been 
applied. 

The results of these calculations are shown in Fig. [TJ The following statements can be made 
concerning the status: 

• M 0u (A) varies slowly with A. This might be due to the fact that only neighboring neutrons 
in a nucleus contribute to the decay [2H E3] • 
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Figure 1. Nuclear matrix element 

calculations for for light neutrino 

exchange. ISM = Interacting shell model |24l 
125] , SRQRPA = self-consistent renormalized 
quasi-particle random phase approximation 
|26j , pnQRPA = proton- neutron quasi particle 
random phase approximation [27], GCM = 
generating coordinate method [28J, IBM = 
interacting boson model [291 EQl [31] (matrix 
elements are scaled by 1.18 to estimate the 
effect if the UCOM short range correlation 
instead of the Jastrow type would have been 
used [9] , in [30] an error of 30% is estimated) , 
pHBF= projected Hartree-Fock-Bogoliubov 
model 132]. 



Figure 2. Ratio of expected 0^/3/3 events per 
kg target mass for the different matrix element 
models normalized to 130 Te. 



• For the ISM, all values are systematically lower than for other methods. Possible reasons 
for this effect are discussed in the literature |244 [S3] ■ 

• For a given isotope the calculations spread by typically a factor of 2-3, i.e. a factor of 4-9 
for Ty 2 . Some groups estimate a possible range for M 0u (see indicated errors in Fig. jlj). 
However it is unclear what to quote as confidence interval for theoretical calculations. 

• The role of short range correlations has been studied and the UCOM correction has emerged 
as favorable [35j. Alternatively, a self consistent implementation was first applied to 
SRQRPA [33] and later to other methods [361 32| and resulted in small changes. 

• Experimental input can cause a sizable shift of the result. For example charge exchange 
reaction measurements of 150 Nd( 3 He,t) and 150 Sm(t, 3 He) |37] result in a quenching factor 
of 0.75 for the gA coupling and hence a reduction of the matrix element by 25% for 150 Nd 
|38j . In this calculation, deformation was treated for the first time in a QRPA calculation. 
For 76 Ge and 76 Se, the proton and neutron valence orbital occupancies have been measured 
[Ml EI]- If the models are adjusted to reproduce these values, the ISM result increases by 
15% [36J while the QRPA results are reduced by about 20% |41| I42j. Hence the difference 



between ISM and QRPA becomes half as large. 



In order to see whether some isotopes are better suited for Ov/3/3 searches from a theoretical 
point of view, the number of expected decays for a given exposure (in units of kg-yr) can be 
compared using Eqs. [8] and 10 without the factors e and r]. The A dependent parameters are 
then the phase space factor, the matrix element and the molar mass. For the comparison it is 
sufficient to look at the ratio of decay rates and in this case, some of the systematic effects of the 
matrix element calculations cancel since there are typically correlations among the isotopes for a 
given method, [j Fig. [2] shows these ratios for the different models normalized to the decay rate 
of 130 Te. One sees that 76 Ge is less favorable. The expected decays per kg vary between 20% and 
50% of the rate of 130 Te. In other words: if all experimental parameters (number of background 
events, efficiency, etc.) were the same then one would need a factor of 2-5 more target mass in 
a 76 Ge experiment to have the same sensitivity. In reality, all experimental parameters like the 
energy resolution and the background have to considered as well. 



5. Current experimental situation 

In the last decade mainly three experiments contributed to OvfiP searches. Heidelberg-Moscow 
was for more than a decade the most sensitive one and reported evidence for this decay. The 
two others were Cuoricino and NEMO-3. Before, IGEX reported with an exposure of 8.8 kg-yr 
a limit of 1®J 2 > 1.57 • 10 25 yr for 76 Ge (90% C.L.)[43j. Currently EXO-200, Kamland-Zen, and 
GERDA are taking data. All experiments are discussed in this section. 



5.1. The Heidelberg- Mo scow experiment 

The Heidelberg- Moscow experiment operated between 1990 and 2003 five germanium detectors 
made out of isotopically enriched material (~86% 76 Ge, 11 kg). The diodes were mounted in 
copper cryostats with copper, lead, and polyethylene shielding. The total exposure was 71.7 kg-yr 
and the average count rate in the interval 2-2.1 MeV was about 0.17 cnts/(keV-kg-yr) (for the 
period 1995-2003). The energy resolution (full width at half maximum, FWHM) was about 
3.5 keV at Qpp which is the best value of all Qvj3j3 experiments. Part of the collaboration finds 
evidence for a peak at Qpp with 28.75 ±6.86 events which converts to Ty 2 = (1-19^0 23) 

• 10 25 yr 

[6]. Note that only a statistical error is quoted. Another study finds that e.g. extending the 
energy window used in the data fit increases this background and hence decreases the signal 
count by up to 40% (Tab. 3.8 and 4.6 of reference [44]). 

In a later publication [JJ the claim was strengthened by a pulse shape analysis which 
preferentially selected Ov(3(3 events due to their localized energy deposition in the detectors. 
Backgrounds from gammas with multiple Compton scatterings exhibit different pulse shapes. 
The background is reduced to a surprisingly low level of ~0.015 cnts/(keV-kg-yr). The final fit 
reports a yield of 11.32±1.75 signal events and the ratio 11.32/1.75 = 6.5 is called the significance 
of the peak (Fig. 9b in [7]). The signal yield was then converted to Ty 2 = (2.23±g;|f) • 10 25 yr. 
There are several problems with this analysis. 

• The fit error on the signal count is too small. The smallest 68% Poisson credibility interval 
is between 8.1 and 15.2 for a probability distribution which peaks at 11.3, i.e. a factor of 2 
larger than the quoted interval. Due to the existing (small) background the ±1<t interval 
should become even larger. 

• The probability that the background (~ 2.2 events in the central 3 keV of the peak) 
fluctuates to the observed number of 13 events or more is 5T0~ 7 which converts to a 

2 This can be seen for example from Tab. 1 of reference [26]. For every isotope, 20 different calculations are listed 
which vary by typically a factor of 2 for a given A. If 20 ratios M 0v (A) / M 0v (A') for two isotopes A and A' are 
calculated, the variation is reduced to 30%. 



significance of about 5<r. Systematic effects like the uncertainty of the background might 
reduce this value. 

• In the conversion to T®J 2 using Eq. (|8j), an efficiency e of 100% is used although no value is 
explicitly quoted. All but three events in the peak are part of an earlier selection (labelled 
"HNR+NN" in [7]). For the latter the efficiency was 62% |45j . Hence one expects also for 
this analysis a value much smaller than 100%. 

The central T®J 2 value and the errors are consequently not correct in [7] and the significance is 
smaller than quoted although still high. 

5.2. The Cuoricino experiment 

Cuoricino [16] operated 62 Te02 crystals with a total mass of 40.7 kg (11.3 kg of 130 Te) between 
2003 and 2008. At a temperature of < 10 mK, the heat capacitance is very low and an energy 
deposition inside a crystal results in an increase of typically 0.1 mK/MeV which is measured with 
Neutron Transmutation Doped germanium thermistors. The latter have a resistance of 100 MO 
and show a strong temperature dependence which converts to a dependence of 3 Mfi/MeV. 
This bolometric technique has been proven to work with a good resolution FWHM of typically 
6-10 keV at 2.6 MeV. 

The total 130 Te exposure was 19.75 kg-yr. The background at Qpp was 0.17 cnts/(keV-kg-yr) 
if normalized to the total mass, i.e. similar to Heidelberg-Moscow. No signal was found and a 
lower limit of T°J 2 > 2.8 • 10 24 yr (at 90% C.L.) was set for 130 Te. This limit is not sensitive 
enough to scrutinize the Heidelberg-Moscow result. 

5.3. The NEMOS experiment 

In NEMO-3 |47] thin foils made out of 7 different /3/3 isotopes (9 kg in total) were located in 
a drift chamber with a magnetic field. Outside of the drift region was a calorimeter made out 
of plastic scintillator blocks with photo multiplier tube (PMT) readout (FWHM for electrons 
15% / y/ E[MeV]) . The 0vf3f3 reconstruction efficiency is only about 8%. On the other hand, 
the topological event reconstruction largely reduces backgrounds from locations other than the 
source foil as well as internal decays with gammas or alphas. Only 2vf3f3 events with poor energy 
reconstruction can not be discriminated. The background is about 1.2 • 10 -3 cnts/(keVkg-yr) 
at 3 MeV (« Qpp of 100 Mo and 82 Se). For all isotopes 2i>/3/3 half lives are reported with 
impressive signal to background ratios of up to 76. For 0^/3/3 decay, the 90% C.L. limits are 
Tfj 2 > 1.0 • 10 24 yr for 100 Mo and Ty 2 > 3.2 • 10 23 yr for 82 Se. Again, NEMO-3 is not sensitive 
enough to confirm or reject the Heidelberg-Moscow claim. 

5.4- The Kamland-Zen experiment 

In Kamland-Zen [3] a balloon of 1.54 m radius made out of 25 //m thick nylon is inserted into the 
Kamland detector and filled with xenon doped scintillator (~290 kg 136 Xe). The energy (FWHM 
~10% at Qpp) and position of the decay (resolution a ~ 15 cm/ y/E(MeV)) is reconstructed with 
PMTs located at a radius of 9 m which cover 34% of the solid angle. 

Since the scintillator is very pure and can be doped easily with several % xenon, this 
experiment has the largest target mass and lowest background - if normalized to the total 
massj^] However the energy resolution is the poorest. 

Kamland-Zen started data taking in 2011. An unexpected background peak at about 2.6 MeV 
is dominating the spectrum around Qpp = 2.458 MeV. It limits the experimental sensitivity and 
the best explanation for the origin is 110m Ag (7i/ 2 = 250 d) cosmogenically produced in 136 Xe 
while it was above ground. This background is expected to be reduced to a negligible level 

3 The Xe mass fraction is ~2.5% and the fiducial volume is ~43%. 



by an ongoing scintillator purification campaign. After an exposure of 78 days first results are 
T% > 5.7 • 10 24 yr (at 90% C.L.) and Tffe = (2.38 ± 0.14) • 10 21 yr [3]. Like for Cuoricino and 
NEMO-3 this limit is not yet sensitive enough to scrutinize Heidelberg-Moscow. 

Additional 700 kg of Xe with 90% enrichment are available by the end of 2012 and are 
expected to be deployed in a cleaner balloon in the near future. Kamland-Zen is therefore 
expected to be the first experiment with ton scale isotope mass. In the more distant future it is 
foreseen to improve the energy reconstruction with light collectors attached to the PMTs and a 
new liquid scintillator [4*8] . 



5.5. The EXO-200 experiment 

EXO-200 |2j operates a liquid xenon TPC of 40 cm diameter and 40 cm length (175 kg of liquid 
Xe, 100 kg fiducial mass). A wire plane in the middle is biased to -8 kV such that electrons 
drift to one of the two ends where the (x, y) position of the electron cloud is reconstructed with 
2 planes of wires and the total charge is measured. The produced scintillation light is detected 
with large area avalanche photo diodes behind the charge collecting wires. 

Since there is a strong anti-correlation between the light amplitude and the ionization signal, 
combining both into a new quantity greatly improves the energy resolution to FWHM ~ 3.9% at 
Qpp. The background is at a quite low level of ~ 1.5 • 10 -3 cnts/(keV-kg-yr) after a geometrical 
cut on the electron cloud is applied to discriminate Compton scattered photon events. 

EXO has recently published results on the 2v(3(3 half life Tffc = (2.1 ± 0.2) • 10 21 yr and a 

limit on the OvfiP half life 1% > 1.6 • 10 25 yr (at 90% C.L.) of 136 Xe 0. 

This result can be used to test the claim of Heidelberg-Moscow. For a given matrix element 
calculation the expected number of signal events for EXO can be estimated from Eq. [8 and 10 
The experimental numbers of EXO (32.5 kg-yr exposure, 55% reconstruction efficiency j 80.6% 
enrichment) and Heidelberg-Moscow (28.75±6.86 events, 71.7 kg-yr exposure, efficiency 100%, 
enrichment 86%) enter. EXO-200 reports ra b s = 1 (^obs = 5) events in an energy window of 
±1(7 (±2(7). 

A Bayesian formulation can be applied with a null hypothesis (H: EXO observes only 
background) and an alternative hypothesis (H: EXO observes a 0uf3f3 signal with scaled 
Heidelberg- Moscow event counts). The priors for the expected number of signal events tt s 
and for background events iTb are Gaussian and the priors for both hypothesis are set to 
n(H) = tt{H) = 0.5. Mean and sigma of ir s are listed in Tab. [2] Mean and sigma of the 
background are 4.1±0.3 (7.5 ± 0.7) events in the ±1(7 (±2<t) energy window. 

The posterior probability is then 

n(m= n(H)-p(D\H) 

P{ ' ir(H) ■ p(D\H) + n(H)p(D\H) 1 ' 



with 



p{D\H) = J ir s {x)ir b (y)P{n ohs \x + y)dxdy (12) 
p(D\H) = [ 7r b (y\h^b)P(n bs\y)dy (13) 



Here -P(n b s |y) is the Poisson function for mean y and n b s observed events. The results for 
p{H) are listed in Tab. [5j For the Tiibingen-Bratislava calculations (labelled "QRPA max" and 
"QRPA min") the maximum and minimum of 20 ratios M 0l/ ( 136 Xe)/M 0l/ ( 76 Ge) are calculated 

An efficiency is not explicitly given in [5] but an effective value can be estimated using Eq.[8]and the information 
from the publication. 



Table 2. Bayesian posterior probabilities p(H) using EXO-200 data for the hypothesis that 
the Ois/3/3 signal of Heidelberg- Moscow is correct. Probabilities are given for different matrix 
element calculations and for the ±1<t and ±2cr energy windows. 

method in ±1 a window in ±2 a window 



expected signal evts. p(H) in % expected signal evts. p(H) in % 



QRPA max 


4.4 ± 1.1 


4 


6.1 ± 1.5 


6 


QRPA min 


2.8 ±0.7 


11 


3.9 ±0.9 


16 


ISM 


10.6 ±2.5 


0.1 


14.8 ±3.5 


0.2 


GCM 


14.3 ±3.4 


0.03 


19.9 ±4.8 


0.05 


pnQRPA 


6.3 ±1.5 


1 


8.8 ±2.1 


2 


IBM 


6.1 ± 1.5 


1 


8.6 ±2.1 


2 



from the values of Tab. 1 of reference [26J. This procedure was suggested by one of the authors 
|49j . For the other models the ratios of the central values shown in Fig. [l]are taken. None of 
the results strengthen the hypothesis of a Oz^/3/3 signal. The ISM and GCM calculations strongly 
disfavor them while the exclusion for the Tubingen-Bratislava calculations are not very strong. 
The reported EXO-200 data were taken within a period of 7 months. Hence more stringent 
statements are expected soon. 

EXO-200 is approved to run for 4 more years. For a following phase the spectroscopic 
identification of the daughter nucleus of Xe is foreseen. In this case the experiment will be 
background free. Only 2vf3j3 events with poorly reconstructed energy can obscure a signal. 

5.6. The GERDA experiment 

GERDA uses the germanium detectors of Heidelberg-Moscow and IGEX and - in a second phase 
- new ones. The f3f3 emitter mass is about 13 kg in the first phase. The detectors are supported 
by a minimal amount of material with low radioactivity in a 4 m diameter cryostat filled with 
liquid argon. Argon serves as cooling medium and shield against external radioactivity. The 
latter is complemented by 3 m of water which is instrumented with PMTs to veto background 
from muons by the detection of their Cherenkov light. 

GERDA started commissioning in 2010 and found an unexpected large background from 42 Ar 
which could be reduced by avoiding electrical fields around the detectors and by an encapsulation 
of the diodes. Since November 2011 the first phase of data taking is ongoing. The background 
is at the level of 0.02 cnts/(keV-kg-yr) and hence almost an order of magnitude smaller than the 
equivalent number of Heidelberg-Moscow. Due to a data blinding procedure no result on 0^/3/3 
is expected before spring 2013. A preliminary result for 2uf3f3 of Ty 2 = (1.88 ± 0.10) • 10 21 yr is 
reported [50] . 

A second phase will start early 2013 with additional new detectors (~18 kg of 
76 Ge). The background is expected to be reduced from currently 0.02 cnts/(keV-kg-yr) to 
0.001 cnts/(keV-kg-yr) due to a liquid argon instrumentation and a different detector type with 
enhanced pulse shape discrimination power. 

6. Other future experiments 

The past and running experiments have been discussed above but there are others under 
construction or as R&D efforts. They use additional isotopes, and various other detection 
mechanisms and background reduction methods, see Tab. [3j 

One important experimental parameter is the fraction rj of the 0v(3(3 isotope (see Eq. ^ and 
Tab.fll. 130 Te is the only one with a large natural abundance. For all other elements enrichment 



of the OvPfl isotope is mandatory. 136 Xe is easiest to enrich with gas centrifuges since this is 
the heaviest isotope and it is already a gas. Other materials like 76 Ge can be converted to a gas 
(GeF4) and then processed. For calcium and neodymium this path is currently not available 
and R&D on alternative methods is ongoing. 

Another important number is Qpp. Larger values are not only better because of larger 
G 0l/ but also because the background from natural decay chains falls off fast beyond 2.6 MeV. 
Experiments using 76 Ge have to compensate by careful material selection and good energy 
resolution. 

The experiments can be grouped into 2 classes. In calorimetric experiments only the total 
energy (ionization or scintillation) is measured. In tracking experiments the two electrons are 
measured independently, i.e. the angular distribution between the electrons is also known. The 
latter is interesting to study the origin of the underlying physics in case 0vf3f3 is observed. 

In the following, the experiments under construction are discussed. 

CUORE is a continuation of Cuoricino with close to 1000 Te02 crystals of 750 g each (in 
total ~ 200 kg of 130 Te). The crystal production is almost finished and all major hardware items 
are ready or close to. A first tower with 52 crystals (CUORE-0) has been assembled and data 
taking in the Cuoricino cryostat was scheduled for July 2012 |51j . 

After a commissioning phase CUORE is expected to start in 2015. 

Majorana will operate similar germanium diodes like GERDA (~ 27 kg of 76 Ge) in vacuum 
in a compact cryostat made out of electro-formed copper. This self-made copper is expected 
to have a factor > 100 smaller thorium, uranium and radium contaminations compared to 
commercial copper such that the background index is about 0.001 cnts/(keV-kg-yr). The 
shielding is completed by commercial copper, lead, and polyethylene. Operations should start 
in 2013. 

In a later phase a combined GERDA and Majorana germanium experiment with order ton 
scale mass is envisioned. 

CANDLES operates 96 scintillating CaF 2 crystals (0.3 kg of 48 Ca) in a liquid scintillator. 
Both are in an acrylic container inside a water tank with PMTs. Data taking started in 2011 
but the sensitivity is limited due to the small target mass. Successful R&D on 8 Ca enrichment 
is crucial for this approach. 

NEXT is a high pressure xenon gas TPC (1.1 m diameter, 1.4 m length, pressure 10-15 
bar, mass 90-130 kg 136 Xe). The time of the 0v(3(3 decay is determined by the detection 
of the (primary) scintillation light with PMTs. The deposited energy and event topology 
is reconstructed from the ionization signal. The drifting electron cloud passes at the end a 
volume of higher electric field such that the electrons are moderately accelerated. Consequently, 
they can excite xenon but not ionize it. The resulting (secondary) scintillation light (electro- 
luminescence) is proportional to the number of electrons and detected with the same PMTs (for 
energy reconstruction) and with a plane of SiPMs (for position reconstruction). The energy 
resolution FWHM is expected to be < 1% at Qpp and the event topology will allow to reject 
backgrounds very effectively to a level of 8 • 10 cnts/(keV-kg-yr). 

Construction will start in 2013 and physics data taking in 2015. 

SuperNEMO is a planned continuation of NEMO-3 with much improved performance 
(factor 4 reconstruction efficiency, factor 2 in energy resolution, factor 6 in background). A 
demonstrator module with 7 kg of 82 Se is under construction and expected to start data taking 
in 2014. 

SNO+ dissolves wl ton nat Nd in 780 tons of liquid scintillator (44 kg of 150 Nd). The 
scintillator is in a 12 m diameter acrylic vessel which is surrounded by pure water in an 
18 m diameter water tank. The latter holds the 9500 8 inch PMTs for light detection. An 
energy resolution FWHM of ~7% at Qpp is expected with an extremely low background of 
< 10~ 6 cnts/(keV-kg-yr) if normalized to the total scintillator mass. The main background 



Table 3. Select 


ion of experiments. 






experiment 


isotope 


mass [kg] 


method 


start / end 


ref. 




past experiments 






Heidelberg-Moscow 


76 Ge 


11 


ionization 


-9003 




Cnoricino 


130 Te 


11 


bolometer 


-2008 


146 


NEMO-3 


100 Mo, 82 Se 


7,1 


track. +calorim. 


-2011 


1471 




current experiments 






EXO-200 


136 Xe 


175 


liquid TPC 


2011- 


121 


f\ m 1 a n rl - ^ pti 


136 Xe 


330 


liquid scintil. 


2011- 


131 


GERDA-I/ GERDA-II 


76 Ge 


15/35 


ionization 


901 1-/ 901 3- 


50 


CANDLES 


48 Ca 


0.35 


scint. crystal 


901 1- 


59 




funded experiments 






NEXT 


136 Xe 


100 


gas TPC 


901 5 


53 


CuoreO/ Cuore 


i30 Te 


10/200 


bolometer 


901 1-1 901 5- 


54 


Majorana Demo. 


76 Ge 


30 


ionization 


2013 


[55J 


SuperNEMO demo. /total 


82 Se 


7/100 


track. +calorim. 


2014-/?? 


[56 


SNO+ 


i50 Nd 


44 


liquid scint. 


2013 


m\ 




proposal, proto-typing 






Cobra 


116 Cd 




solid TPC 




[58J 


Lucifer 


82 Se 




bolom. +scint. 




[59j 


DCBA/MTD 


i50 Nd 


32 


tracking 




[60j 


MOON 


82 Se, 100 Mo 


30-480 


track. +scint. 




m\ 


AMoRE 


100 Mo 


100 


bolom. +scint. 




[62j 


Cd exp. 


116 Cd 




scint. 




[63^ 



is 2vf3f3 due to the poor energy resolution and the relatively short half life T 2 ^ of 150 Nd (see 
Tab. [l]). Scintillator filling and doping with Nd is expected for 2013. Since 150 Nd can not be 
enriched effectively at the moment, i] is small which limits the SNO+ sensitivity. 

Lucifer is an R&D effort which investigates scintillating crystals with (3/3 emitters like ZnSe 
which are operated as bolometers. The simultaneous detection of phonons and photons allows 
to identify backgrounds from e.g. surface events which are expected to dominate in CUORE. 

Other ongoing R&D efforts are not discussed here. References are listed in Tab. [3j 

7. Comparison of experiments 

For a comparison of the sensitivities of the experiments a relative scaling factor for the 
different matrix elements and phase spaces has to be applied. This factor can be estimated 
from Fig. [§) The values used here are f A (Ge) = 0.35(0.2 - 0.5), f A (Se) = 1.1(0.7 - 1.4), 
f A (Mo) = 2.1(1.1 - 3.2), f A (Te) = 1, / A (Xe) = 0.55(0.4 - 0.7) and / A (Nd) = 1.2(0.4 - 2.0). 
The numbers in parentheses are the full range. 

If the number of background events is large, Eq. ^ can be used to estimate the experimental 
sensitivity. A relative figure-of-merit can then be defined as 

FOM = f A .e. V . ] l^ E (14) 

The relative T®J 2 sensitivity scales with the live time t of an experiment like FOM-\/t. Tab.jijiists 
the performance numbers of the experiments discussed above. "Kamland-Zen2" is the improved 
experiment after the purification of the scintillator (assumed factor 5 smaller total background) 



Table 4. Comparison of relative figure-of-merit (FOM), lower half life limit Ty 2 after 4 yr 
live time, and resulting upper limit on mpp. For mpp, the entire range of matrix element values 
including the indicated error bars in Fig. [T] are used, /a is the average scale factor for a given 
isotope taken from Fig. [2] AE is the energy window which is taken to be 1(2) FWHM for 
experiments with > 0.5% (< 0.5%) resolution. Note that the efficiency is reduced by 0.7 if 
AE = 1-FHWM. FOM is defined in the text. Masses are total masses or fiducial masses. The 
background and enrichment fraction has to be scaled accordingly. 



exp. 


mass 
Heel 


fA 


bkg. 

r 10~ 3 cnt l 
LkcV-kg-yr-l 


AE 
[keV] 


eff. 


enrich. 


FOM 


1/ 1 

10 25 vr 


mpp 








past experiments 










Hd-Moscow 


11 


0.35 


120 


7 


1 


0.86 


1 


1.9 


170-530 


Cuoricino 


41 


1 


170 


16 


0.9 


0.28 


1 


0.4 


210-500 


NEMO-3 


6.9 


2.1 


1.2 


400 


0.06 


0.9 


0.3 


0.1 


310-900 








running 


; experiments 










EXO-200 


100 


0.55 


1.5 


100 


0.55 


0.81 


6 


4.2 


75-170 


Kami. -Zen 


12800 


0.55 


0.05 


250 


0.31 


0.023 


4 


2.6 


90-220 


Kaml.-Zen2 


12800 


0.55 


0.01 


250 


0.31 


0.06 


22 


15 


40-90 


GERDA-I 


15 


0.35 


20 


8 


0.8 


0.86 


2 


3.9 


120-370 


GERDA-II 


35 


0.35 


1 


6 


0.85 


0.88 


20 


18 


60-170 








experiments under construction 








Major. -Dem. 


30 


0.35 


1 


6 


0.9 


0.9 


20 


17 


60-170 


CUORE 


750 


1 


10 


12 


0.9 


0.27 


19 


7.5 


50-110 


SNO+ 


780000 


1.5 


0.0002 


230 


0.33 


5.6E-5 


3 


0.8 


100-240 


NEXT 


100 


0.55 


0.8 


25 


0.25 


0.9 


9 


5.2 


70-160 








proposed experiments 










S.NEMO 


100 


1.1 


0.1 


200 


0.2 


0.9 


14 


6.9 


55-140 


Lucifer 


100 


1.1 


1 


10 


0.9 


0.5 


50 


19 


33-85 



and the upgrade to one ton xenon mass. For comparison, the FOM numbers, the expected 90% 
C.L. limits for 4 yr of live time, and the corresponding mpp limits are given. For the latter, 
the entire spread of the matrix elements of Fig. [T] including the error bars are used. 

For running (and past) experiments the achieved performance values are used which might 
improve with time. For the others the anticipated performance numbers are taken. 

As a graphical representation, the relative sensitivity of the experiments as a function of live 
time is shown in Fig. [3] This value is calculated from Eq. Q by 



fou fA-e-V-M-t 

H / 2> *(B-AE-M-? [ ' 



Here \I/(Abkg) is the "average" 90% C.L. upper limit of the number of signal events for Abkg 
background events calculated according to the method discussed in [M] . 

A few comments should be made concerning the interpretation of Tab. [4] and Fig. [3j 

• The factor Ja and hence FOM has a full spread (not a) of ~ rfc 30%-70%. Thus the curves 
in Fig. [3] could be replaced by bands which would make the figure however unreadable. 

• The sensitivities discussed here are calculated for Oz//3/3 exclusion limits. For a positive signal 
claim, the situation is different. A good energy resolution like the ones for germanium or 
bolometer experiments will allow to identify a narrow line at the correct energy. This is 
extremely valuable if the existence of the rarest ever observed decay will be claimed. 




Figure 3. Relative experimental sensitivity for the 0^/3/3 half life limit versus live time for 
different experiments. 



• To estimate relative sensitivities for mpp, the inverse square root has to be taken of the 
curves shown in Fig. [3] and the FOM numbers in Tab [4} The variation due to the spread of 
the matrix elements is reduced in this case as can be seen from the last column of Tab. |4j 
i.e. the spread for the lower or upper value of the mpp interval is reduced. 

If instead limits for other particle physics parameters are calculated, the scaling will be 
different. For heavy neutrino exchange, the mass limit scales with y ^1/2 ■ The factor Ja 
will change in this case but for most models the variations is less than 30% |651 [66] . 

• The 2vf3p background is irreducible and can only be avoided with an energy resolution 
a < 1 — 2% at Qpp. This requirement depends of course strongly on Ty 2 which varies by 
a factor of 300 for the isotopes considered. 

• Of the ongoing experiments, Kamland-Zen(2) should have the largest potential. However 
a low background bolometer experiment like Lucifer or an improved CUORE experiment 
will be even more sensitive. The sensitivity of them would still grow almost linearly after 
4 yr (see Fig. [3j> . 

• Germanium experiments can be competitive to e.g. current xenon experiments with a factor 
of 3 more mass despite the fact that the phase space factor is small. Required is however a 
factor of 10 lower background than the one of GERDA-I. 

• Systematic effects like the uncertainty on the fiducial volume or the knowledge of the 
background level are not taken into account here. 

The goal of many searches is to reach a sensitivity equivalent to mag ~ 19 meV which would 



cover practically the entire expected range for the inverted neutrino mass hierarchy. For 76 Ge, 
this corresponds to half lives of (1.5 — 15) • 10 27 years if the entire span of matrix elements of 
Fig. [T] is taken into account. These values should be compared to the expected sensitivity of 
GERDA-II or Majorana Demonstrator of about 1.7-10 26 y. This demonstrates that exploring the 
entire mass band of the inverted hierarchy is a long term enterprise^] The equivalent numbers 
for 136 Xe, 82 Se, and 130 Te are (0.6 - 3.5) • 10 27 yr, (0.6 - 3.8) • 10 27 yr, and (0.5 - 2.5) • 10 27 yr, 
respectively. They seem to be easier to reach if one compares them to the T®J 2 limits in Tab. |4j 

Considering the example of the Heidelberg-Moscow claim and the uncertainties of the matrix 
elements, one can conclude that a signal has to be observed in several isotopes to establish Oz^/3/3. 

To improve the credibility of the result a blind analysis should be performed. This technique 
is nowadays standard in particle physics experiments and should be adopted in this field as well. 

It is also worth mentioning that continuously new ideas for 0vf3f3 experiments are coming up 
|67j and the one who will make a discovery might not be listed in this publication. 



8. Summary 

Neutrinoless double beta decay violates lepton number and the experimental programs are 
therefore on equal footing to proton decay searches. It might also be the only practical process 
which allows to test whether neutrinos are Majorana particles. The motivation for several large 
efforts in this field is therefore obvious. 

For a long time, the Heidelberg-Moscow experiment has dominated the field and its claim of 
a 0vj3j3 signal has not been scrutinized since 2001. The recent EXO-200 limit does not support 
this claim but can not refute it due to the spread of the nuclear matrix element calculations. 
GERDA does not suffer from such uncertainties and will unblind the data in spring 2013. Then 
the combined data from GERDA, EXO-200 and Kamland-Zen should be sensitive enough for a 
meaningful test. 

Beyond this next step, experiments want to reach a sensitivity to explore the mpp region 
of the inverted neutrino mass hierarchy. This will eventually require ton scale experiments. 
Kamland-Zen will be the first one but whether the entire range will be covered will depend on 
the achievable background level. 

For a convincing claim of a 0v/3/3 signal a good energy resolution is important and the 
detection with several isotopes. 
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